Linearized f(R) gravity: Gravitational radiation and Solar System tests 
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We investigate the linearized form of metric /(_R)-gravity, assuming that f(R) is analytic about 
R — so it may be expanded as f(R) = R + a-zR 2 /2 + . . . . Gravitational radiation is modified, 
admitting an extra mode of oscillation, that of the Ricci scalar. We derive an effective energy- 
momentum tensor for the radiation. We also present weak-field metrics for simple sources. These 
are distinct from the equivalent Kerr (or Schwarzschild) forms. We apply the metrics to tests that 
could constrain f(R). We show that light deflection experiments cannot distinguish /(i?)-gravity 
from general relativity as both have an effective post-Newtonian parameter 7 = 1. We find that 
planetary precession rates are enhanced relative to general relativity; from the orbit of Mercury 
we derive the bound |<22 < 1.2 x 10 18 m 2 . Gravitational-wave astronomy may be more useful: 
considering the phase of a gravitational waveform we estimate deviations from general relativity 
could be measurable for an extreme-mass-ratio inspiral about a 10 6 Mq black hole if 1 0,2 1 > 10 17 m 2 , 
assuming that the weak-field metric of the black hole coincides with that of a point mass. However 
Eot-Wash experiments provide the strictest bound \a?\ <2x 1CP 9 m 2 . Although the astronomical 
bounds are weaker, they are still of interest in the case that the effective form of f(R) is modified 
in different regions, perhaps through the chameleon mechanism. Assuming the laboratory bound is 
universal, we conclude that the propagating Ricci scalar mode cannot be excited by astrophysical 
sources. 

PACS numbers: 04.50.Kd, 04.25.Nx, 04.30.-w, 04.70.-s 



I. INTRODUCTION 

General relativity (GR) is a well tested theory of grav- 
ity [l[; so far no evidence has been found that suggests 
it is not the correct classical theory of gravitation. How- 
ever, there are many unanswered questions that remain 
regarding gravity which motivate the exploration of alter- 
nate theories: What are the true natures of dark matter 
and dark energy? How should we formulate a quantiz- 
ablc theory of gravity? What drove inflation in the early 
Universe? Is GR the only theory that is consistent with 
current observations? Moreover, the majority of the tests 
that have been carried out to date have been in the weak- 
field, low-energy regime ELGl : in the laboratory 0, 3 , 
within the Solar System [5|, IfJ or using binary pulsars jj . 
It is not unreasonable to suppose that GR would begin 
to break down at higher energies. 

Over the coming decade, a new avenue for testing rela- 
tivity will be opened up, through the detection of gravita- 
tional waves (GWs) using the existing ground-based GW 
detectors, the Laser Interferometer Gravitational- Wave 
Observatory (LIGO) Virgo and GEO [HG1, 

and the proposed space-based GW detector, the Laser 
Interferometer Space Antenna (LISA) [l3|, These 
detectors will observe GWs generated during the inspi- 
ral and merger of binary systems comprising one or more 
black holes (BHs). The GWs are generated in the strong- 
field regime, while the components are highly relativis- 
tic and the spacetime is evolving dynamically: GW as- 
tronomy will open a new window into the strong-field 
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regime of gravity, c omp lementing traditional electromag- 
netic observations [I3. A comparison of the GWs ob- 
served from such systems with the predictions of GR will 
provide powerful tests of the theory in a region yet to be 
explored. 

The radiation generated during the final merger and 
ringdown of two BHs will offer tests of GR in the high- 
est energy and most dynamical sector, but it is thought 
that the most sensitive tests will come from LISA ob- 
servations of extreme-mass-ratio inspirals (EMRIs) [l6| . 
An EMRI involves the inspiral of a stellar-mass compact 
object, a white dwarf, neutron star or BH, into a massive 
BH in the centre of a galaxy. The mass of the compact 
object is typically 1-10M Q , while the mass of the mas- 
sive BH (for sources in the LISA band) will be ~ 10 5 - 
10 7 M Q , so the mass-ratio is of the order of ~ 10~ 7 -10~ 4 . 
This extreme mass-ratio means that the inspiral proceeds 
slowly, and on short timescales the compact object acts 
like a test particle moving in the background spacetime 
of the central BH. LISA will detect ~ 10 5 cycles of grav- 
itational radiation generated while the compact object is 
in the strong field of the spacetime, and this encodes a 
detailed map of the spacetime structure outside the cen- 
tral BH. This idea was first elucidated by Ryan [13, [U| , 
who showed that, for an arbitrary stationary and axisym- 
metric spacetime in GR, the multipole moments of the 
spacetime enter at different orders in an expansion of the 
frequency of small vertical or radial oscillations of circu- 
lar, equatorial orbits. As these frequencies are in princi- 
ple observable in the GWs generated during an inspiral, 
it should be possible to measure the multipole moments 
from an EMRI observation and hence test whether the 
central object is a Kerr BH: according to the no-hair 
theorem, a Kerr BH is described completely by its mass 
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M and spin angular momentum J jl9H23f . and its mass 
multipole M; (M = M) and mass-current multipole mo- 
ments Si (Si = J) are determined from these according 

to m 



Mi +iSi=M 



J 

M 



(1) 



The multipole expansion is not a convenient way to char- 
acterize arbitrary spacetimes, since the Kerr metric itself 
requires an infinite number of multipoles to fully char- 
acterize. Subsequent authors have instead adopted the 
approach of considering bumpy BH spacetimes (25l - |28| . 
which deviate from the Kerr metric by a small amount 
and depend on some parameter, e, such that e = is pre- 
cisely the Kerr solution. Relatively small perturbations 
to the Kerr solution can be detected in EMRI observa- 
tions due to small differences in the precession frequen- 
cies that accumulate over the 100 000 waveform cycles 
that will be detected. There are also certain qualita- 
tive features that could be smoking-guns for a departure 
from the Kerr metric, such as ergodicity in the orbits j28| . 
persistent resonances [29[ or a shift in the frequency of 
plunge |1 Hp. 

The majority of the work to date has focused on space- 
times that are solutions in GR, but which deviate from 
the Kerr solution. However, if GR was not the correct 
theory of gravity, this could also lead to detectable signa- 
tures in the observed gravitational waves. Certain alter- 
native theories of gravity, including f(R), do admit the 
Kerr metric as a solution, since it has vanishing Ricci ten- 
sor, = dHHH. However, the Kerr metric need not 
be the expected end state of gravitational collapse [33| . If 
a Kerr BH existed in an alternative theory, the geodesies 
would be the same, but the energy flux carried by the 
GWs could still be different, and so differences would 
show up in the rate of inspiral; although in many cases 
these differences do not appear at leading order. In most 
cases, however, either the Kerr metric is not admitted 
as a solution, or it is not the correct metric to describe 
collapsed objects [HJ . Waveform differences then show 
up as a result of the differences in the instantaneously- 
geodesic orbits of the compact object involved in the 
EMRI. Since the leading-order energy-momentum tensor 
of the GWs often takes the same form as in GR [34[ , this 
is the primary effect and means the problem of testing 
alternative theories through EMRI observations is equiv- 
alent to the spacetime mapping programme within GR 
described previously. 

As a consequence of the difficulties of solving for GW 
emission in alternative theories, work on testing alter- 
native theories of gravity using LISA EMRIs has so far 
been restricted to a few cases. In Brans-Dicke gravity, in 
which the gravitational field is coupled to a scalar field, 
differences show up due to a modification to the inspi- 
ral rate that arises from dipole radiation of the scalar 
field (HJ. Neutron star EMRIs are required since the 
dipole radiation depends on a sensitivity difference be- 
tween the two objects, and the sensitivity is the same for 



all BHs. Lower mass central BHs provide the most pow- 
erful constraints, but a LISA observation of a neutron 
star EMRI into a 10 4 Af Q BH could place constraints on 
the Brans-Dicke coupling parameter that are competi- 
tive with Solar System constraints [35J. In dynamical 
Chern-Simons modified gravity, the action is modified 
by a parity- violating correction, inspired by string the- 
ory [36|, |33|. In this case, the BH solution differs from 
the Kerr solution at the fourth multipole, I = 4, but 
the energy-momentum tensor of gravitational radiation 
takes the same form as in GR [38|]. LISA observations 
of EMRIs should place constraints on the Chern-Simons 
coupling parameter that are an order of magnitude bet- 
ter than will be possible from binary pulsar observations, 
although a full analysis accounting for parameter degen- 
eracies has not yet been carried out [381 ] - 

In this work, we focus our attention on metric f(R)- 
gravity, in which the Einstcin-Hilbert action is modified 
by replacing the Ricci scalar R with an arbitrary func- 
tion f(R). This is one of the simplest extensions to stan- 
dard GR [39|, It has attracted significant interest 
because the flexibility in defining the function f(R) al- 
lows a wide range of cosmological phenomena to be de- 
scribed Hi], S3]. F° r example, Starobinsky Q suggested 
that a quadratic addition to the field equations could 
drive exponential expansion of the early Universe [3]: 
inflation in modern terminology. In this model f(R) = 
R — R 2 / (6T 2 ); the size of the quadratic correction can be 
tightly constrained by considering the spectrum of cur- 
vature perturbations generated during inflation |4(| . 
Using the results of the Wilkinson Microwave Anisotropy 
Probe J!?], HH] , the inverse length-scale can be constrained 
to T ~ 3 x 10" 6 (50/A r )Zp 1 HE Hi, where N is the num- 
ber of e-folds during inflation and lp is the Planck length. 

We consider simple f(R) corrections within the frame- 
work of linearized gravity, and explore what constraints 
LISA might be able to place on the form of f(R) (we will 
not consider cosmological implications where terms be- 
yond linear order could play a significant role). We will 
see that, although the field equations for /(JO-gravity do 
admit the Kerr metric as a solution [3l], |33| , this is not 
necessarily the metric that describes the exterior of col- 
lapsed objects. We consider the modifications to geodesic 
orbits in the weak- field of the f(R) spacetime exterior 
to massive objects and, assuming this also describes the 
weak-field external to a BH, we estimate how observable 
the differences in the precession frequencies will be by 
LISA. We will also describe Solar System and laboratory 
constraints that can be placed on the same model. The 
overall conclusion is that LISA could place constraints 
on /(i?)-gravity, which may be more powerful than those 
in the Solar System, but not as powerful as constraints 
from laboratory experiments. However, the LISA ob- 
servations will probe a different energy scale, so these 
constraints will still be important, particularly if we re- 
gard f(R) as an effective theory that could be different 
in different regimes. 

This paper is organised as follows. We begin with a 
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review of the f(R) field equations. In Sec. Mil we derive 
the linearized equations and in Sec. HVI we apply these to 
find wave solutions. These results can be used to study 
how gravitational radiation is modified for /(i?)-gravity. 
They are largely known in the literature, but are worked 
out here ab initio; they are included as a compendium 
of useful results within a consistent system of notation. 
To be able to accurately model gravitational waveforms 
one needs to know how an object will inspiral. Accord- 
ingly, we derive an effective energy-momentum tensor 
for gravitational radiation in Sec. [V] following the short- 
wavelength approximation of Isaacson [H HJ ■ In Sec. IVII 
we look at the effects of introducing a source term and de- 
rive the weak-field metrics for a point source, a slowly ro- 
tating point source, and a uniform density sphere, recov- 
ering some results known for quadratic theories of gravity. 
These are used in Sec. IVIII to compute the frequencies of 
radial and vertical epicyclic oscillations about circular- 
equatorial orbits in the weak-field, slow-rotation metric, 
and hence to construct an estimate of the dctectability 
of the f(R) deviations in LISA EMM observations. For 
comparison, in Sec. IVIII1 we describe the constraints on 
/(i?)-gravity that can be obtained from Solar System and 
laboratory tests. We conclude in Sec. IIXI with a summary 
of our findings. 

Throughout this work we will use the timelike sign 
convention of Landau and Lifshitz [52j : 

1. The metric has signature (+, — , — , — ). 

2. The Riemann tensor is defined as R^vao = 



vp 



3. The Ricci tensor is defined as the contraction 

Rp,u — R p\v • 

Greek indices are used to represent spacetime indices [i = 
{0, 1, 2, 3} (or \x = {t, r, 8, </>}) and lowercase Latin indices 
are used for spatial indices i = {1,2,3}. Natural units 
with c = 1 will be used throughout, but factors of G will 
be retained. 



II. DESCRIPTION OF /(_R)-GRAVITY 

A. The action and field equations 

General relativity may be derived from the Einstcin- 
Hilbert action [53, l53fl 



Seh[<7] 
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(2) 



In f(R) theory we make a simple modification of the 
action to include an arbitrary function of the Ricci scalar 
R such that [53 



S[g] = 



i 



16ttG 



f(R)^9d 4 x. 



(3) 



Including the function f(R) gives extra freedom in defin- 
ing the behaviour of gravity. While this action may 
not encode the true theory of gravity it could contain 
sufficient information to act as an effective field theory, 
correctly describing phenomenological behaviour (55| ; it 
may be that as an effective field theory, a particular 
f(R) will have a limited region of applicability and will 
not be universal. We will assume that f(R) is analytic 
about R = so that it can be expressed as a power se- 
ries 0, [II [1111 



f{R) = a Q 



aiR+a i R2 + a i R " 



(4) 



Since the dimensions of f{R) must be the same as of R, 
[a n ] = [R] (1 - n \ To link to GR we will set ai = 1; any 
rescaling can be absorbed into the definition of G. 

Various models of cosmological interest may be ex- 
pressed in such a form, for example, the model of 
Starobinsky ^ 



f(R) = R + XR 



can be expanded as 



f(R) = R 



An 

R 



■R 2 



R 2 

1 + W 2 



Xn(n + 1) 
2*8 



R 4 



(5) 



(6) 



Consequently such a series expansion can constrain 
model parameters, although we cannot specify the full 
functional form from only a few terms. 

The field equations are obtained by a variational prin- 
ciple; there are several ways of achieving this. To derive 
the Einstein field equations from the Einstein-Hilbert ac- 
tion one may use the standard metric variation or the 
Palatini variation [53[. Both approaches can be used for 
f(R), however they yield different results [H, S3- 
lowing the metric formalism, one varies the action with 
respect to the metric g^ u . Following the Palatini formal- 
ism one varies the action with respect to both the metric 
and the connection r p Mt/ , which are treated as inde- 
pendent quantities: the connection is not the Levi-Civita 
metric connection^ 

Finally, there is a third version of /(i?)-gravity: mctric- 
affine /(iJ)-gravity ■ This goes beyond the Pala- 

tini formalism by supposing that the matter action is de- 
pendent on the variational independent connection. Par- 
allel transport and the covariant derivative are divorced 
from the metric. This theory has its attractions: it allows 



1 Imposing the condition that that the metric and Palatini for- 
malisms produce the same field equations, assuming an action 
that only depends on the metric and Riemann tensor, results in 
Lovelock gravity Lovelock gravities require the field equa- 

tions to be divergence free and no more than second order; in 
four dimensions the only possible Lovelock gravity is GR with a 
potentially nonzero cosmological constant I60H62H . 
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for a natural introduction of torsion. However, it is not 
a metric theory of gravity and so cannot satisfy all the 
postulates of the Einstein equivalence principle [l[ : a free 
particle does not necessarily follow a geodesic and so the 
effects of gravity might not be locally removed [59| . The 
implications of this have not been fully explored, but for 
this reason we will not consider the theory further. 

We will restrict our attention to metric /(i?)-gravity. 
This is preferred as the Palatini formalism has undesir- 
able properties: static spherically symmetric objects de- 
scribed by a polytropic equation of state are subject to 
a curvature singularity [io|, HH, Ull • Varying the action 
with respect to the metric g^" produces 



(7) 



where □ = g^V^Vt, is the d'Alembertian and a prime 
denotes differentiation with respect to R. Proceeding 
from here requires certain assumptions regarding sur- 
face terms. In the case of the Einstein-Hilbcrt action 
these gather into a total derivative. It is possible to sub- 
tract this from the action to obtain a well-defmcd varia- 
tional quantity foTi Ir38j . This is not the case for general 
f(R) [69(- However, since the action includes higher- 
order derivatives of the metric we are at liberty to fix 
more degrees of freedom at the boundary, in so doing 
eliminating the importance of the surface terms (39l . [70l | . 
Setting the variation SR. = on the boundary allows us 
to subtract a term similar to in GR [7l|. Thus we have 
a well-defined variational quantity, from which we can 
obtain the field equations. 

The vacuum field equations are 

f'Rfiv - V M V„/' + g^nf - L 9ltv = 0. (8) 
Taking the trace of our field equations gives 



f'R + 3\3f - 2/ = 0. 



(9) 



If we consider a uniform flat spacetime R = 0, this equa- 
tion gives [EH 



a = 0. 

In analogy to the Einstein tensor, we define 
/ / / / 

so that in a vacuum 
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0. 



(10) 

(11) 

(12) 



Acting upon this with the covariant derivative 

i 



+ /'V [R^ - -R 9fiv 
- (DV. - V„D) /'. 



(14) 



The second term contains the covariant derivative of the 
Einstein tensor and so is zero. The final term can be 
shown to be 

(□V„ - V„D) /' = R„u V/', (15) 

which is a useful geometric identity (72[. Using this 

8nGV^T^ = 0. (16) 

Consequently energy-momentum is a conserved quantity 
in the same way as in GR, as is expected from the sym- 
metries of the action. 

III. LINEARIZED THEORY 

We start our investigation of f(R) by looking at lin- 
earized theory. This is a weak-field approximation that 
assumes only small deviations from a flat background, 
greatly simplifying the field equations. Just as in GR, 
the linearized framework provides a natural way to study 
GWs. We will see that the linearized field equations will 
reduce down to flat-space wave equations: GWs are as 
much a part of /(i?)-gravity as of GR. 

Consider a perturbation of the metric from flat 
Minkowski space such that 



(17) 



where, more formally, we mean that h^ v = eH^ v for a 
small parameter e0 We will consider terms only to 0(e). 
Thus, the inverse metric is 



(18) 



B. Conservation of energy- momentum 



where we have used the Minkowski metric to raise the 
indices on the right, defining 

= rTv vp fhr P . (19) 

Similarly, the trace h is given by 

h = rf v h„ v . (20) 

All quantities denoted by "h" are strictly 0{e). 
The linearized connection is 

r (1) V = \if x {d^h Xv + d v h Xll - a A v)« ( 21 ) 



If we introduce matter with a stress-energy tensor T^ v , 
the field equations become 



Guv = 



(13) 



2 It is because we wish to perturb about flat spacetime that we 
have required f(R) to be analytic about R = 0. 



5 



To O(e) the covariant derivative of any perturbed quan- 
tity will be the same as the partial derivative. The Rie- 
mann tensor is 



In GR we use the trace-reversed form where = 
— {hj2)r]^ v . This will not suffice here, but let us look 
for a similar solution 



(22) 

where we have raised the index on the differential opera- 
tor with the background Minkowski metric. Contracting 
gives the Ricci tensor 

R {1 \v = \{d^d P K + dvdph' - d^duh - Uhfiu), (23) 

where the d'Alembertian operator is □ = if^d^d^. Con- 
tracting this with rf v gives the first-order Ricci scalar 



R {1) = d^dphw - Dh. 



(24) 



To O(e) we can write f(R) as a Maclaurin series 

f(R) = a Q + R {1) - 1 (25) 
f(R) = l + a 2 R il) . (26) 

As we are perturbing from a Minkowski background 
where the Ricci scalar vanishes, we use (TTU)) to set ao = 0. 
Inserting these into pT|) and retaining terms to 0(e) 
yields 

Q (1 \u = R {1 \,~d^(a 2 R^)+r,^D(a2R (1) )-^V^- 

(27) 

Now consider the linearized trace equation, from (J9j) 

gW = R^+3D(a 2 R^)~2R^ 

G (1) = 3D(a 2 .R (1) )-i? (1) , (28) 

where Q^ 1 ' = rf^Q^ 1 p. v - This is the massive inhomoge- 
neous Klein-Gordon equation. Setting Q — 0, as for a 
vacuum, we obtain the standard Klein-Gordon equation 

+T 2 i? (1) = 0, (29) 

defining the reciprocal length (squared) 



1 

3a 2 



(30) 



For a physically meaningful solution T 2 > 0: we con- 
strain f(R) such that a 2 < [z3-[zH| - From T we define 
a reduced Compton wavelength 



A R = - 



(31) 



associated with this scalar mode. 

The next step is to substitute in h pu to try to find 
wave solutions. We want a quantity h^y that will satisfy 
a wave equation, related to h^ v by 



(32) 



(33) 



The only rank- two tensors in our theory are: h^v, rj^ v , 
R^'ptvi and dpd v ; h^ u has been used already, and we wish 
to eliminate R^ 1 '^, so we will try the simpler option 
based around r]^. We want B^ v to be 0(e); since we 
have already used h, we will try the other scalar quantity 
Therefore, we construct an ansatz 



ba 2 R {1) 



(34) 



where a 2 has been included to ensure dimensional con- 
sistency and b is a dimensionlcss number. Contracting 
with the background metric yields 

h = Aba 2 R {1) - h, (35) 

so we can eliminate h in our definition of h pv to give 



ba 2 R {l) - ~ ) rip,. 



Just as in GR, we have the freedom to perform a gauge 
transformation IZZl: the field equations are gauge- 
invariant since we started with a function of the gauge- 
invariant Ricci scalar. We will assume a Lorenz, or de 
Dondcr, gauge choice 



V/i M „ = 0; 



or for a flat spacetime 



d»Ku = 0. 



(37) 
(38) 



Subject to this, from (|23|) . the Ricci tensor is 



1 1 flV 



2bd^(a 2 R^) + D[ V- -r hl 



+ £ (jR (l) + S(l) )V/ 



(39) 



Using this with ((25J in ((27]) gives 



0« 



- W Id I 



G 



-(6 + 1) 



/ ' ' ■-' 



d^d^R^ + ^R^v^ 



(40) 



Picking b = — 1 the final term vanishes, thus we set [76j, 



a 2 



RW 



h 



h 



hpv = - ( a 2 i? (1) + - ] 



(41a) 
(41b) 
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From (l24l) the Ricci scalar is 



in this case 



= □ ^r(i) _ ^ _ □(-4a 2J R( 1 ) - 71) 

= 3D(a 2 i? (1) ) + hjfi. (42) 
For consistency with (j2"5)l , we require 

- = (43) 

Inserting this into (|40|) . with 6 = — 1, we see 

- lav = S ( V; (44) 

we have our wave equation. 

Should a 2 be sufficiently small that it can be regarded 
an 0(e) quantity, we recover the usual GR formulae to 
leading order within our analysis. 



IV. GRAVITATIONAL RADIATION 

Having established two wave equations, (f28|) and (|l4"]l . 
we now investigate their solutions. Consider waves in a 
vacuum, such that Q^ v = 0. Using a standard Fourier 
decomposition 



hfiv — hfwikp) exp yikpX^) , 
= R(q p ) exp (iq p xP) , 



(45) 
(46) 



where fc M and are four-wavevectors. From (|44[) wc 
know that k^ is a null vector, so for a wave travelling 
along the z-axis 

fc M =w(l, 0,0,1), (47) 
where w is the angular frequency. Similarly, from (1281) 

g" = f n, 0, 0, V^ 2 - T 2 ) , (48) 



for frequency fl. These waves do not travel at c, but have 
a group velocity 



yn 2 - t 2 
n ' 



(49) 



provided that T 2 > 0, v < 1 = c. For < T, we find an 
evanescently decaying wave. The travelling wave is dis- 
persive. For waves made up of a range of frequency com- 
ponents there will be a time delay between the arrival of 
the high-frequency and low-frequency constituents. This 
may make it difficult to reconstruct a waveform, should 
the scalar mode be observed with a GW detector (79| . 

From the gauge condition (|38|) we find that fc M is or- 
thogonal to h^, 



k huif — 0, 



(50) 



h 3v = 0. 



(51) 



Let us consider the implications of (|43[) using equations 
M and dSU, 



□ ( 4a 2 R {1) 



-h 
Oh 







(52) 



For nonzero RW (as required for the Ricci mode) there 
is no way to make a gauge choice such that the trace 
h will vanish [Zfl [z^. This is distinct from in GR. It 
is possible, however, to make a gauge choice such that 
the trace h will vanish. Consider a gauge transformation 
generated by £ M which satisfies D£ M = 0, and so has a 
Fourier decomposition 



£ At exp(ifc p x p ) 



(53) 



A transformation 
h 



+ dptu + d v ^ - v^dPZp, (54) 



would ensure both conditions ([38JI and (|44|) are satis- 
fied 1531. Under such a transformation 



(55) 



We may therefore impose four further constraints (one 
for each £ M ) upon /i M „. We take these to be 



h 0v = 0, 



0. 



(56) 



This might appear to be five constraints, however we have 
already imposed (|51[) . and so setting /loo = automati- 
cally implies /103 = 0. In this gauge we have 

V = Ku - a 2 R il) r )l , v , h= -4o 2 i? (1) . (57) 

Thus behaves just as its GR counterpart; we can 
define 



"0 








0" 





h + 


h x 










-h+ 














0. 



(58) 



where h + and h x are constants representing the ampli- 
tudes of the two transverse polarizations of gravitational 
radiation. 

It is important that our solutions reduce to those of GR 
in the event that f(R) = R. In our linearized approach 
this corresponds to a 2 — > 0, T 2 — > oo. We see from (|4"5)l 
that in this limit it would take an infinite frequency to 
excite a propagating Ricci mode, and evanescent waves 
would decay away infinitely quickly. Therefore there 
would be no detectable Ricci modes and we would only 
observe the two polarizations found in GR. Additionally, 
h^ v would simplify to its usual trace-reversed form. 
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V. ENERGY-MOMENTUM TENSOR 

We expect gravitational radiation to carry energy- 
momentum. Unfortunately, it is difficult to define a 
proper energy-momentum tensor for a gravitational field: 
as a consequence of the equivalence principle it is possi- 
ble to transform to a freely falling frame, eliminating the 
gravitational field and any associated energy density at 
a given point, although we can still define curvature in 
the neighbourhood of this point [H, [77}. We will do 
nothing revolutionary here, but will follow the approach 
of Isaacson [50|, lEH- The full field equations (jBJ) have 
no energy-momentum tensor for the gravitational field 
on the right-hand side. However, by expanding beyond 
the linear terms we can find a suitable effective energy- 
momentum tensor for GWs. Expanding Q pv in orders of 
h,. 



'in' 



(2) 



(59) 



We use (B) for the background term instead of (0) to 
avoid potential confusion regarding its order in e. So far 
we have assumed that our background is flat; however, 
we can imagine that should the gravitational radiation 
carry energy-momentum then this would act as a source 
of curvature for the background [|oj]. This is a second- 
order effect that may be encoded, to accuracy of C(e 2 ), 
as 



(60) 



By shifting Q^ 2 \ u to the right-hand side we create an 
effective energy-momentum tensor. As in GR we will av- 
erage over several wavelengths, assu ming that the back- 
ground curvature is on a larger scale }34 l53j. 



'(B) 



Q (2 \ 



(61) 



By averaging we probe the curvature in a macroscopic re- 
gion about a given point in spacctimc, yielding a gauge- 
invariant measure of the gravitational field strength. The 
averaging can be thought of as smoothing out the rapidly 
varying ripples of the radiation, leaving only the coarse- 
grained component that acts as a source for the back- 
ground curvature The effective energy-momentum ten- 
sor for the radiation is 



1 



2(2) 



(62) 



Having made this provisional identification, we must 
set about carefully evaluating the various terms in (|59|). 
We begin as in Sec. IHII bv defining a total metric 



9nu = lnv + h 



in/ ■ 



(63) 



3 By averaging we do not try to localise the energy of a wave to 
within a wavelength; for the massive Ricci scalar mode we always 
consider scales greater than \r. 



where 7 M „ is the background metric. This changes our 
definition for h pv : instead of being the total perturba- 
tion from flat Minkowski, it is the dynamical part of the 
metric with which we associate radiative effects. Since 
we know that Q^ B \ V is 0(s 2 ), we decompose our back- 
ground metric as 



(64) 



where is 0(e 2 ) to ensure that pvp is also 0(e 2 ). 

Therefore its introduction will make no difference to the 
linearized theory 

We will consider terms only to 0(e 2 ). We identify 
T^ P pi/ from (HJ). There is one small subtlety: whether 
we use the background metric or jf v to raise indices 
of <9 M and h pv . Fortunately, to the accuracy considered 
here, it does not make a difference; however, we will con- 
sider the indices to be changed with 7^. We will not 
distinguish between dp and V^ B ^, the covariant deriva- 
tive for the background metric: to the order of accuracy 
required covariant derivatives would commute and V' B ^ 
behaves just like <9 M . Thus 



V 

2 1 



dp (h\ v - a 2 R { 1} 7a 
d v (hxp - a 2 R {1) ixp 
dx (/i M „ - a 2 i? (1) 7, t „ 



(65) 



and 



r (2)p M „ = -^h px (d^h Xu + d v h Xlx - d x v) 



(h pX -a 2 R^h pX ) [d,(h x 



a 2 



i? (1) 7A, 



+ d v yh Xfl - a 2 i? (1) 7A M 
- d\ \h pv - a 2 i? (1) 7^ 



(66) 



For the Ricci tensor we can use our linearized expres- 
sion, (|39)) . for the first-order term, 



per 



(67) 



= a 2 d^R {1) + \r (1) 1^- 

The next term is 

p(3) _ 9r (2)" _ 5r (2) p +y{i) p r (i) c 

x jj,a x pv 



- I -d^hapdji 9 + h ap 



Jpu - 9 u d p ( h 



a 2 J 



-R (1) 7ctm ) _ dpd p [h av - a 2 R 



' ap 



'lav 



d p \i v (d p h afl - d a h pfJ ) - a 2 d° R (1} d a h pv 
a\ (2R {l) 8pd v R {l) + 3dpR {1) d u R (1) 
R^D^R^J)}. (68) 
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The d'Alembertian is = j^d^d u . 

To find the Ricci scalar we contract the Ricci tensor 
with the full metric. To 0(e 2 ), 



and 



R (B) = ^ R (B)^ 

R (2) = r »R(2)^ _ h ^ R W^ 



(69) 



(70) 


f (B) 


= a 2 RW 






= 1 


2a 2 h P "d p d„R {1 ' ) 




= a 2 RW 




= a 2 RW 



(75) 
(76) 
(77) 
(78) 



Using these 



f(B) = R (B) 

f W = R (i) 

/( 2) = R (2) + ?1rM\ 



(71) 

(72) 
(73) 
(74) 



We list a zeroth-order term for clarity. i?( B ' is 0{e 2 ). 
Combining all of these 



= R {2 \, + a 2 ( 7 „D( B ) - d,d v ) R^ IrW^ + °± ( 7 ^ D (B) _ Riif _ 

- a 2lp X P d a d p R (1) + ^d p R w fahpv + 8 v h pp - d p h pv ) + a 2 (rWrW mv + ±R^\^ 



~l P v 



d p R w d n R {1 ' 



(79) 



It is simplest to split this up for the purposes of averag- 
ing. Since we average over all directions at each point, 
gradients average to zero [HI, [77} 

(d P V) = 0. (80) 

As a corollary of this we have 

(Ud p V) = - (Vd p U) . (81) 

Repeated application of this, together with our gauge 
condition, and wave equations, ([2"8]l and (|4"4")l , allows 
us to eliminate many terms. Those that do not average 
to zero are the last three terms in (|T9"|) . plus 



RW 



- A d,K p d v h pa + ^d p R^d v R^ 



R (2)) = (^ R (lf) 



(82) 
(83) 



[rVR^^) = Ur^O^R^ + \l^R {1)2 ) ■ (84) 



Combining terms gives 



j(2) 



fiu ] — \ ^o p n ap o y n ^ u 



Thus we obtain the result 
1 



32ttG 



dnh ap dX° + 6a 2 2 d p R {1) d„R {1 



(85) 



(86) 



In the limit of a 2 — > we obtain the familiar GR result as 
required. The GR result is also recovered if R^ = 0, as 
would be the case if the Ricci mode was not excited; for 
example, if the frequency was below the cutoff frequency 
T. Rewriting the effective energy- momentum tensor in 
terms of metric perturbation h pvi using (1571) . 



t, 



1 



32ttG 



1 



d p h ap d v h f " J + -d p hd„h 



(87) 



These results do not depend upon 03 or higher-order co- 
efficients [34 1 . 

The effective energy-momentum tensor could be used 
to constrain the parameter a 2 through observations of 
the energy and momentum carried away by GWs. Of 
particular interest would be a system with a frequency 
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that evolved from u> < T to to > T, as then we would 
witness the switching on of the propagating Ricci mode. 
If we could accurately identify the cutoff frequency we 
could accurately measure a-i- However, see Sec. IVIII CI 
for further discussion of why this is unlikely to happen. 



We will again assume a Minkowski background, consider- 
ing terms to 0(e) only. To solve the wave equations (|28|) 
and (l44l) with this source term we use a Green's function 



(a + r 2 )& r (x,x') = s(x-x'), 



(89) 



VI. f(R) WITH A SOURCE 



where □ acts on x. The Green's function is familiar as 
the Klein-Gordon propagator (up to a factor of — i) [8l| 



Having considered radiation in a vacuum, we now add 
a source term. We want a first-order perturbation, so the 
linearized field equations are 

Q^\ v = %^GT, V . 



-(x, x) 



d A p exp [— ip ■ (x — x') 
(2tt) 4 T 2 -p 1 



(90) 



This can be evaluated by a suitable contour integral to 



& r (x,x') 



f dw . , ,,-,1 



da; 

— exp [-iuj{t-t)\ 



4irr 
1 

Airr 



exp 



exp 



i [oj — T I r 



T 2 -uj 



cu 2 >T 2 

uj 2 <r 2 



(91) 



where we have introduced t = x , t' = x and r 
la: - x'\. For T = 



%(x,x') = 



S(t -t'-r) 



4nr 



(92) 



the familiar retarded-time Green's function. We can use 
this to solve (gj) 



Ku{x) = -WttG J d 4 x'%{x,x')T^(x') 



(93) 



This is exactly as in GR, so we can use standard results. 
Solving for the scalar mode 

R {1 \x) = SttGT 2 f d 4 x'& T (x,x')T(x'). (94) 



To proceed further we must know the form of the trace 
T(x'). In general the form of i?^^(x) will be complicated. 



A. The Newtonian limit 

Let us consider the limiting case of a Newtonian source, 
such that 

T 00 = p; |Too| » \T 0i \; \T 00 \ > |T^|, (95) 
with a mass distribution of a stationary point source 

p = MS{x'). (96) 
This source does not produce any radiation. As in GR 
4GM 



hoi — hij — 0. 



(97) 



Solving for the Ricci scalar 



(98) 



Combining these in (|41b[) yields a metric perturbation 
with nonzero elements 



hj/i 



2GM 

r 

2GM 



exp(— T?*) 



1 - 



exp(— T?') 



(99) 



Sij. 



Thus, to first order, the metric for a point mass in f{R)- 
gravity is (5(| HI, HU 



ds 2 



2GM 



1 



exp(— T?') 



2GM 



3 

exp(— Tr) 



dt 2 



dS 2 , (100) 



using dE 2 = dx 2 + dy 2 + dz 2 . This is not the linearized 
limit of the Schwarzschild metric (although it is recovered 
as 0,2 —> 0, T —> oo) 84|. This metric has already been 
derived for the case of quadratic gravity, which includes 
terms like R 2 and R^R^ in the Lagrangian [H [73, H, 
l86j |. In linearized theory our f(R) reduces to quadratic 
theory, as to first order /(i?) = R + a 2 R 2 /2. 

Extending this result to a slowly rotating source with 
angular momentum J, we then have the additional 
term [l3 



2G 



Xk, 



(101) 
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where e^ k is the Levi-Civita alternating tensor. The met- 
ric is 

2GM 



ds z 



1 



AG J 



1 



1 



exp(— Tr) 



(xdy — ydx) di 



2GM 



1 



exp(— Tr) 



d£ 2 , (102) 



where z is the rotation axis. This is not the first-order 
limit of the Kerr metric (aside from as ai — > 0, T — > oo). 

In /(i?Wravity Birkhoff's theorem no longer ap- 
plies [85l - l89j : the metric about a spherically symmet- 
ric mass does not correspond to the equivalent of the 
Schwarzschild solution. The distribution of matter in- 
fluences how the Ricci scalar decays, and consequently 
Gauss' theorem is not applicable. Repeating our anal- 
ysis for a (nonrotating) sphere of uniform density and 
radius L 



AGM 



too 



0. 



(103) 



as in GR, and for the point mass, but 



-6GM 



exp(— Tr) 



TLcosh(TL) -sinh(TL) 



TL 3 



-6GM 6XP( Tr) T 2 S(TL), 



(104) 



defining S(TL) in the last lineH The metric perturbation 
thus has nonzero first-order elements 86, 88, 89] 



^oo = - 2GM [1 + exp(-Tr)S(TL)] ; 
hij = - 2GM [1 - cxp(-Tr)S(TL)] Si 



(105) 



where we have assumed that r > L at all stages@ 

Solving the full field equations to find the exact metric 
in f(R) is difficult because of the higher-order derivatives 
that enter the equations. However, we expect a solution 
to have the appropriate limiting form as given above. 

It has been suggested that since R = is a valid so- 
lution to the vacuum equations, the BH solutions of GR 
should also be the BH solutions in f(R) [HI, HI]. How- 
ever, while the Kerr solutions are solutions of the vac- 
uum field equations, the presence of a source complicates 
the issue; it may be that the end point of gravitational 
collapse is not the Kerr solution, and so astrophysical 
BHs in /(i?)-gravity may not be the same as their GR 
equivalents. We have seen that having a nonzero stress- 
energy tensor at the origin, because of (|28p . forces R to 
be nonzero in the surrounding vacuum, although it will 
decay to zero at infinity [9(j. While one cannot gener- 
alise straightforwardly from our simple (5-function sources 



to complete BH solutions, because of the horizon in the 
BH spacetime, these solutions suggest that astrophysical 
BHs could be different from the Kerr solution!! If astro- 
physical BHs are not described by the Kerr metric, these 
weak-field metrics provide a reasonable candidate for the 
alternative form. 

If the astrophysical BHs in /(i?)-gravity have a differ- 
ent structure from their GR counterparts it should be 
possible to distinguish between theories by observing the 
BHs that form. It is this possibility that we focus on in 
the next section. Even in the event that the BH space- 
times do coincide, we could still detect differences in the 
properties of extended sources. 



B. The weak- field metric 



It is useful to transform the weak-field metric, (|102l) . 
to the more familiar form 

AG J 

ds 2 = A(7)dt 2 + — ^- sin 2 dd<pdt - B(7)d7 2 - F 2 df2 2 . 

(106) 

The coordinate r is a circumferential measure, as in the 
Schwarzschild metric, as opposed to r, used in preceding 
sections, which is a radial distance (an isotropic coordi- 
nate) [5J, [9(| • To simplify the algebra we introduce the 
Schwarzschild radius 



r s = 2GM. 



(107) 



In the linearized regime, we require that the new radial 
coordinate satisfies 



1 



r 



1 



exp(— Tr) 



rs 



1 - 



exp(— Tr) 



(108) 
(109) 



This can be used as an implicit definition of r in terms 
of r . To first order in rs jr [90| 



A(r) 



1 



rs 
r 



1 



exp(— Tr) 



(110) 



We see that the functional form of 500 is almost un- 
changed upon substituting r for r; however r is still in 
the exponential. 

To find B(r) we consider, using (11091) . 



— = d In r 

r 



1 + Trgr exp(— Tr)/6r 1 dr 

1 + (r s /2?) [1 - exp(-Tr)/3] J V' 

(111) 



4 H(0) = 1/3 is the minimum of r.(TL). 

5 Inside the source ijW = -(6GM/L 3 )[1 - (TL + 1) cxp(-TL) x 
sinh(Tr)/Tr], 



There is currently no proof of the uniqueness of the Kerr solutions 
as the end state of gravitational collapse in f(R), although there 
does exist a similar result for the closely related Brans-Dicke 
theory HUii. 
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Thus 



form 



df 2 = -7T 



1 + Tr$r exp(— Tr)/6r 
r 2 { 1 + (r s /2r) [1 - exp(-Tr)/3] 



dr 2 . (112) 



dty 
dr 



(117) 



The term in braces is [-B(r)] . We assume that in the 
weak-field 



where w M is the 4- velocity. For a circular equatorial orbit, 
setting (i = r gives the frequency of the orbit wo = d(f>/dt 
as 



e ~ — 
r 



(113) 



is small. Then the metric perturbations from Minkowski 
arc small. Expanding to first order [9^ 



w 



GJ 1 2A' [2GJ\ 



(118) 



in which a dash denotes d/df and the +/— sign denotes 
progradc/rctrogradc orbits. The definition of proper 
time gives 



B(r) = l + r 4 
r 



1 - 



exp(— Tr) 



Trg exp(— Tr) 



(114) 



t = A 



4GJw 



r2, .2 



-1/2 



(119) 



In the limit T — > oo, where we recover GR, Air) and 
5(r) tend to their Kerr (Schwarzschild) forms. 

In the following sections we will use these weak-field 
metrics (in both coordinates) with astrophysical and lab- 
oratory tests of gravity to place constraints on f(R). 



We now have both t and (f> = u)$t as functions of r; sub- 
stitution into (|115[) and (|116|) allows us to find the energy 
and angular momentum in terms of r. 

From the Hamiltonian H = g^u^u" we can obtain the 
general equation of motion for massive particles, using 
the substitutions 



VII. EPICYCLIC FREQUENCIES 

One means of probing the nature of a spacetime is 
through observations of orbital motions [28 1. Wc will 
consider the epicyclic motion produced by perturbing a 
circular orbit. There are two epicyclic frequencies asso- 
ciated with any circular-equatorial orbit, characterizing 
perturbations in the radial and vertical directions respec- 
tively [95| . We will start by deriving a general result for 
any metric of the form of (|106j) . and then specialise to our 
f(R) solution. Wc will work in the slow-rotation limit, 
keeping only linear terms in J. 

An orbit in a spacetime described by (|106[) has as con- 
stants of motion: the orbiting particle's rest mass, the en- 
ergy (per unit mass) of the orbit E and the z-component 
of the angular momentum (per unit mass) L z . Using 
an over-dot to denote differentiation with respect to an 
affine parameter, which wc identify as proper time t, 



. E 

t= A 



2GJ_ 



2GJE 



Ar 3 r 2 sin 2 1 



(120) 
(121) 



where we have linearized in J, as appropriate for the 
slow-rotation limit. With these replacements, the general 
timclikc geodesic equation takes the form 



~ 2 r A2 



AB 



4GJEL 2 



ABr 3 
= V(¥,9,E,L Z ). 



L- 



r 2 sin 2 1 



(122) 



To compute the epicyclic frequency we imagine the orbit 
is perturbed by a small amount, while E and L z are un- 
changcdQ For radial perturbations r = r(l + 6), where r 
is the radius of the unperturbed orbit, the orbit under- 
goes small oscillations with frequency 



i 2 Q 2 



1 d 2 V 



a d 



2 dr - 



(123) 



E = At 
L 



2GJ 



sin 



r 2 sin 2 1 



2GJ 



sin 2 6t. 



(115) 
(116) 



We will consider perturbations of circular-equatorial or- 
bits; orbits such that r = r = 9 = and 8 = ir/2. The 
timelike geodesic equation can be written in the covariant 



Small vertical perturbations 8 = ir/2 + S oscillate with 
frequency 



;2q2 _ _ 
1 "vert — 



lB(r) d 2 V 



dd 2 



(124) 



r,e = Tr/2 



We will denote A(r) = A, B(r) = B, A'(r) = A , etc. 
differentiating the potential from (|122j) wc find 



7 It is not possible for the orbit to be perturbed without changing 
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i'ad 



-E 2 



A 



■12 



A 



A 3 B 



a'b' 



2A 2 B A 2 B 2 



B 



B 



AGJEL. 



3 

^3 



A' 



r B \ A 



a'b 



2A 



A 



AB 2AB 

a'b' b' 2 



B 

2# 



B 



2A 2 B 



A 2 B 2 



6GJEL Z / A 



AB 
■a 



2AB 



A Br 4 



A' f, 



B 

(a 2 b 



B 



B 



2B 2 r 



tt2 



B 



AB 



B 3 r 2 



2 I " ABr 2 



2B 



B r 



3 



(125) 

(126) 
(127) 



To simplify (|125[) we have used conditions imposed by 
setting V and dV/dr equal to zero for circular, equatorial 
orbits. These results hold for any metric of the general 
form (|106p . subject to the slow-rotation condition, which 
we have used to linearize in J at various stages. 

A. Gravitational-wave constraints 

We are interested in whether or not the deviation aris- 
ing from the f(R) correction would be observable. In 
principle, the deviations will be observable if the orbit 
looks sufficiently different from orbits in the Kerr metric@ 
To quantify the amount of difference, we need to iden- 
tify orbits between the two spacetimcs, and for circular- 
equatorial orbits there is a natural way to do this: by 
identifying orbits with the same frequency uyk since this 
is a gauge-invariant observable quantity [96|. The quan- 
tity 

A(w ,T) = n(uj ,T) - Q(lo 0i T -> oo) (128) 

characterizes the rate of increase in the phase difference 
between the f(R) trajectory and the Kerr trajectory with 
the same frequency and spin parameter A physical 
effect is in principle observable if it leads to a signif- 
icant phase shift in a gravitational waveform over the 
length of an observation. Thus, a simple criterion for 
the f(R) theory to be distinguishable from GR would 
be that T Q b s A > 27r, for observation period T b s - This 
is a significant oversimplification, since we have assumed 
that only the orbital frequency has been matched to a 



the energy or angular momentum. However, these corrections 
are quadratic in the amplitude of the perturbation, and so we 
can ignore them at linear order. 

Here we assume that the end point of gravitational collapse is 
not the Kerr solution, and that the weak-field f(R) metric is a 
reasonable approximation to the true astrophysical BH solution. 
If it were Kerr, the epicyclic frequencies would not differ between 
f(R) and GR. 

By comparing the trajectory to the T — y oo limit of the trajectory 
rather than the exact Kerr result ensures that we are taking the 
same slow rotation limit in both cases, and will not be dominated 
by 0(J 2 ) corrections. 



Kerr value, while small changes in the other parameters 
such as the BH mass and spin, the orbital eccentricity 
and inclination, and so on, could mimic the effects of an 
f(R) deviation. On the other hand, we are also keeping 
the orbital frequency fixed whereas we will observe inspi- 
rals, and this tends to break the parameter degeneracies. 
Since we are interested in cxtremc-mass-ratio systems, 
for which the inspiral proceeds slowly, it is likely that 
we are being over-optimistic, so these results can be con- 
sidered upper bounds on what could be measurable. A 
fuller analysis accounting for parameter correlations and 
inspiral is beyond the scope of this paper. 

The timcscalc of the systems we are considering is set 
by the BH mass, and the quantities Mujq and MA are 
mass- independent. A duration of a typical EMRI obser- 
vation with LISA will be of the order of a year, and so 
the criterion for detectability becomes 

GMA = 9.8 x 1(T 7 ( f ) [ —\ . (129) 

\ioni Q ) \T ohs J v ; 

In Fig. [T] we show the region of T-wo parameter space 
in which f(R) gravity could be distinguished from GR, 
as defined by this criterion. Each curve represents a 
particular choice for GMA, and the region below the 
curve is detectable in an observation characterized by 
that choice for MA. Equation (|129p indicates that the 
curve GMA = 10~ 6 is what would be achieved in a 
one-year observation for a 10 6 M Q mass BH. The curves 
GMA = 10~ 5 /10 -7 are the corresponding results for a 
1O 7 /1O 5 M mass BH, while the curve GMA = 3x 1CT 7 
represents what would be achieved in a three-year ob- 
servation and so on. We show results for two different 
choices of spin, a = J/(GM 2 ) = and a = 0.5, and 
it is clear that there is not too much difference between 
the two; although the vertical epicyclic frequency is only 
measurable for b / since it coincides with the orbital 
frequency for a = because of the spherical symmetry 
of the potential. The results for the radial epicyclic fre- 
quency do not differ hugely between a = and a = 0.5 
in this weak-field metric approximation. We note also 
that we show results only for prograde orbits. For a^O, 
we can also compute results for retrograde orbits, and 
these differ from the prograde results but only by a small 
amount which is almost indistinguishable on the scale of 
these plots. 
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FIG. 1. Region of parameter space in which f(R) theories can be distinguished from GR when the central BH has spin a = 
(left panel) or a — 0.5 (right panel). Each curve corresponds to a particular specification of the detectability criterion given in 
(|129l) in the text, as identified in the key. Dashed lines correspond to measurements of the vertical epicyclic frequency, while 
solid lines represent measurements of the radial epicyclic frequency. The region below the curve could be distinguishable in a 
LISA observation with that detectability value. 



Our conclusion from Fig. [T]is that, broadly speaking, 
we would be able to distinguish spacetimes with GMT < 
1, for a 1O 6 M BH this corresponds to T < 1CT 9 m _1 . 
Somewhat larger values are accessible at higher frequen- 
cies, but this conclusion must be treated somewhat cau- 
tiously, as the inspiral would pass through that region 
fairly quickly, and those orbits correspond to relatively 
small values of the orbital radius at which the approxima- 
tions that we made deriving the weak-field metric begin 
to break down. For this criterion, the radial epicyclic fre- 
quency is always a more powerful probe than the vertical 
epicyclic frequency. This is to be expected, since the lat- 
ter is generally smaller in magnitude and so fewer cycles 
accumulate over a typical observation. 



VIII. SOLAR SYSTEM AND LABORATORY 
TESTS 

A. Post-Newtonian parameter 7 

The parametrized post-Newtonian (PPN) formalism 
was created to quantify deviations from GR [3, 0. It 
is ideal for Solar System tests. The only parameter we 
need to consider here is 7, which measures the space- 
curvature produced by unit rest mass. The PPN metric 
has components 

<7oo = 1-2*7; 9ij =-(l + 2 7 C/)%, (130) 

where for a point mass 

U(r) = ™. (131) 
r 

The metric must be in isotropic coordinates @, HH . The 
f(R) metric (|100[) is of a similar form, but there is not a 



direct correspondence because of the exponential^ It 
has been suggested that this may be incorporate d by 
changing the definition of the potential U [40l.l57l. IfJOl IfJS 
then 



7 



3 — exp(— Tr) 
3 + cxp(— Tr) ' 



(132) 



As T — > 00, the GR value of 7 = 1 is recovered. How- 
ever, the experimental bounds for 7 are derived assuming 
that it is a constant 0. Since this is not the case, we 
will rcderive the post-Newtonian, or 0(e), corrections to 
photon trajectories for a more general metric. We define 



ds 2 



P(r)dt 2 - Q(r) (dx 2 + dy 2 + dz 2 ) 



(133) 



To post-Newtonian order, this has nonzero connection 
coefficients 



r° 



Qi 



2r ' 
Q'(S^x k 



-pi _ 
1 00 - 



P'x i 
- S jk x l ) 



(134) 



2r 



The photon trajectory is described by the geodesic equa- 
tion 



dx" dx p 



0, 



da 2 ' vp da do 
for affinc parameter a. The time coordinate obeys 



(135) 



da 2 



dx" dx p 
' da da 



0. 



(136) 



10 Our f(R) theory is equivalent to a Brans-Dicke theory with a 
potential and parameter ojbd = [97I [9§L We cannot use the 
familiar result 7 = (1 + ljbd)/(2 + lobd) \M as this was derived 
for Brans-Dicke theory without a potential Q. 
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so wc can rewrite the spatial components of (|135p using 
t as an affine parameter [2[ 



dx u dxP 
~~dT~dT 



dt 2 ' V~ VP " " P dt 
Since the geodesic is null we also have 
dx** dx" 



9nu- 



= 0. 



dt dt 

To post-Newtonian accuracy these become 



(137) 



(138) 



dV 




dx 




dt 2 


[ 2r 2r 


dt 





P' - 



-x 



dx dx' 

"dT"dT' 



= P-Q 



dx 



dt 



(139) 
(140) 



The Newtonian, or zeroth-order, solution of these is prop- 
agation in a straight line at constant speed Q 

x l N =n l t; \n\ = l. (141) 

The post-Newtonian trajectory can be written as 

x l = nH + x pN (142) 

where x l pN is the deviation from the straight line. Sub- 
stituting this into (|139[) and (|140[) gives 

^jf- = -\V{P -Q) + n-V{P- Q)n, (143) 
d:Ep N P - 



— dr"^ < 144 » 

The post-Newtonian deviation only depends upon the 
difference P-Q. From ([TUP)) 

4GM 



P(r) - Q(r) = 



r 

-4U(r) 



(145) 

This is identical to in GR. The result holds not just for 
a point mass, we see, using (|41b[) , 

P(r) — Q(r) = hoo + ha (no summation) 

= h aa + hu, (146) 

and since h^ u obeys (|44|) exactly as in GR, there is no 
difference. We conclude that an appropriate definition 
for the post-Newtonian parameter is 



7 



.900 + Qii 

ZU 



- 1 



(no summation). 



(147) 



Using this, our f(R) solutions ha ve j = 1. This agrees 
with the result found by Clifton [58( 03 Consequently, 
/(i?)-gravity is indistinguishable from GR in this respect 
and is entirely consistent with the current observational 
value of 7 = 1 + (2.1 ± 2.3) x 1CT 5 [H, [|. We must use 
other experiments to put constraints upon f(R). 



11 Clifton [H also gives PPN parameters = 1, £l = 0, £3 = 
and £4 = 0, all identical to in GR. 



B. Planetary precession 

We can also use the epicyclic frequencies derived in 
Sec. lVIII for the classic test of planetary precession in the 
Solar System. Radial motion perturbs the orbit into an 
ellipse. The amplitude of our p erturbation S gives the 
eccentricity e of the ellipse 100]. Unless ljq = f2 ra d the 
epicyclic motion will be asynchronous with the orbital 
motion: there will be precession of the periapsis. In one 
revolution the ellipse will precess about the focus by 



— 



2tt 



ad 



(148) 



where ujq is the frequency of the circular orbit, given in 
(|118[) . The precession is cumulative, so a small devia- 
tion may be measurable over sufficient time. Taking the 
nonrotating limit, the epicyclic frequency is 



"rad 



1 



3r s 



C(T,r s ,r) 



(149) 



defining the function 

exp(— T?*) 



1 



r \ r 



T 2 r 2 cxp(-Tr) 
3 + (l + Tr) exp(-Tr) 
exp(— Tr) 



(150) 



This characterizes the deviation from the Schwarzschild 
case: the change in the precession per orbit relative to 
Schwarzschild is 



An7 = w — zug 



(151) 
(152) 



using the subscript S to denote the Schwarzschild value. 
To obtain the last line we have expanded to lowest order, 
assuming that £ is small0 Since £ > 0, the precession 
rate is enhanced relative to GR. 

Tablc|I]shows the orbital properties of the planets. We 
will use the deviation in perihelion precession rate from 
the GR prediction to constrain the value of £, and hence 
T and 02- All the precession rates are consistent with 
GR predictions (Aw = 0) to within their uncertainties. 
Assuming that these uncertainties constrain the possi- 
ble deviation from GR we can use them as bounds for 
the f(R) corrections. Table IPTl shows the constraints for 
T and 02 obtained by equating the uncertainty in the 
precession rate <7A ra with the f(R) correction, and sim- 
ilarly using twice the uncertainty 2a A ro . The tightest 
constraint is obtained from the orbit of Mercury. Adopt- 
ing a value of T > 5.3 x 10~ 10 m _1 , the cutoff frequency 
for the Ricci mode is > 0.16 s _1 . Therefore it could lie in 



12 There is one term in f that is not explicitly O(e). Numerical 
evaluation shows that this is < 0.6 for the applicable range of 
parameters. 
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TABLE I. Orbital properties of the eight major planets and Pluto. We take the semimajor orbital axis to be the flat-space 
distance r, not the coordinate r. The eccentricity is not used in calculations, but is given to assess the accuracy of neglecting 
terms 0(e 2 ). 





Qpm ! tti a ] r\v ayi q 101 

O -jllllllld 1 Ul ClAlo | 1U1 


Orhital neriod flOll 


Pr'pppQtiion ra1"p 1 1 ? 


P,pppn_~rip_ N Dl 1 


r^iaiiei 


r /i nil 

/ / lu III 


(27r/^o)/yr 


Azu zb <j Aro /m__s yr 


(":■ 


Mercury 


0.57909175 


0.24084445 


-0.040 ± 0.050 


0.20563069 


Venus 


1.0820893 


0.61518257 


0.24 ± 0.33 


0.00677323 


Earth 


1.4959789 


0.99997862 


0.06 ± 0.07 


0.01671022 


Mars 


2.2793664 


1.88071105 


-0.07 ± 0.07 


0.09341233 


Jupiter 


7.7841202 


11.85652502 


0.67 ± 0.93 


0.04839266 


Saturn 


14.267254 


29.42351935 


-0.10 ± 0.15 


0.05415060 


Uranus 


28.709722 


83.74740682 


-38.9 ± 39.0 


0.04716771 


Neptune 


44.982529 


163.7232045 


-44.4 ± 54.0 


0.00858587 


Pluto 


59.063762 


248.0208 


28.4 ±45.1 


0.24880766 



TABLE II. Bounds calculated using uncertainties in planetary perihelion precession rates. T must be greater than or equal to 
the tabulated value, 1 02 1 must be less than or equal to the tabulated value. 

Using (TA ra Using 2a A m 



Planet 


T/10" 11 m" 1 


M/10 18 m 2 


T/10" 11 m" 1 


M/10 18 m 


Mercury 


52.6 


1.2 


51.3 


1.3 


Venus 


25.3 


5.2 


24.6 


5.5 


Earth 


19.1 


9.1 


18.6 


9.6 


Mars 


12.2 


22 


11.9 


24 


Jupiter 


2.96 


380 


2.87 


410 


Saturn 


1.69 


1200 


1.63 


1200 


Uranus 


0.58 


9800 


0.56 


11000 


Neptune 


0.35 


28000 


0.33 


31000 


Pluto 


0.26 


49000 


0.25 


55000 



the upper range of the LISA frequency band (i~3l . [lij or in 
the LIGO /Virgo frequency range |8l-ilG] . The constraints 
are not as tight as those which could be placed using 
gravitational-wave observations. However, as we will see 
in Sec. I VIII C[ it is possible to place stronger constraints 
on T using laboratory experiments. 



Fifth-force tests 



From the metric (|100[) we see that a point mass has a 
Yukawa gravitational potential [U 



V(r) 



GM 



exp(— TV) 



(153) 



Potentials of this form are well studied in fifth-force 
tests [l|, [3, H| which consider a potential defined by a 
coupling constant a and a length-scale A such that 



V(r) 



GM 



1 + a exp 



(154) 



We are able to put strict constraints upon our length- 
scale Ak, and hence 02, since our coupling constant cxr = 
1/3 is relatively large. This can be larger for extended 
sources: comparison with (|105[) shows that for a uniform 
sphere a_R = E(TL) > 1/3. 



The best constraints at short distances come fro m the 
Eot -Wash experiments, which use torsion balances |l03l . 
llO These constrain A_r < 8 x 1CU 5 m. Hence we deter- 
mine |ci2 1 <2x 10~ 9 m 2 . A similar result was obtained 
by Naf and Jetzer [Hj]. This would mean that the cut- 
off frequency for a propagating scalar mode would be 
> 4 x 10 12 s _1 . This is much higher than expected for 
astrophysical objects. 

Fifth- force tests also permit Xn to be large. This de- 
generacy can be broken using other tests; from Sec. IVIII 
we know that the large range for A_r is excluded by plan- 
etary precession rates. This is supported by a result of 
Naf and Jetzer [83| obtained using the results of Gravity 
Probe B @. 

While the laboratory bound on Xr may be strict 
compared to astronomical length-scales, it is still much 
greater than the expected characteristic gravitational 
scale, the Planck length lp. We might expect for a nat- 
ural quantum theory that 02 ~ 0(lp); however lp — 
2.612 x 10 -70 m 2 , thus the bound is still about 60 or- 
ders of magnitude greater than the natural value. The 
only other length-scale that we could introduce would be 
defined by the cosmological constant A. Using the con- 
cordance values A = 1.26 x 1CU 52 m~ 2 ; we see that 
A -1 ^> |a 2 |. It is intriguing that if we combine these two 
length-scales we find Zp/A 1 / 2 = 1.44 x 10~ 9 m 2 , which 
is of the order of the current bound. This is likely to be 
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a coincidence, since there is nothing fundamental about 
the current level of precision. It would be interesting to 
sec if the measurements could be improved to rule out a 
Yukawa interaction around this length-scale. 



IX. SUMMARY AND CONCLUSIONS 

We have examined the possibility of testing f(R) type 
modifications to gravity using future gravitational-wave 
observations and other measurements. We have seen that 
gravitational radiation is modified in /(i?)-gravity as the 
Ricci scalar is no longer constrained to be zero and, in 
linearized theory, there is an additional mode of oscilla- 
tion, that of the Ricci scalar. This is only excited above a 
cutoff frequency, but once a propagated mode is excited, 
it will carry additional energy-momentum away from the 
source. The two transverse GW modes are modified from 
their GR counterparts to include a contribution from the 
Ricci scalar, see (|41a[) , which will allow us to probe the 
curvature of the strong-field regions from which GWs 
originate. However, further study is needed in order to 
understand how the GWs behave in a region with back- 
ground curvature, in particular, when R is nonzero. 

From linearized theory we have deduced the weak- 
field metrics for some simple mass distributions and 
found they are not the BH solutions of GR. Additionally, 
Birkhoff's theorem no longer applies in /(i?)-gravity. If 
the end point of gravitational collapse is not the Kerr so- 
lution, LISA observations of extreme-mass-ratio inspirals 
will be sensitive to small differences in the precession fre- 
quencies of orbits, as small differences lead to secular de- 
phasings that accumulate over the 100 000 waveform cy- 
cles LISA will observe. By computing epicyclic frequen- 
cies for the weak-field, slow-rotation metric we were able 
to estimate the constraints that might come from such 
observations. These indicated that deviations would only 
be detectable when |a-2 1 ^ 10 17 m 2 , assuming an extreme- 
mass-ratio binary with a massive BH of mass ~ 10 6 M©. 
We also discussed constraints that could be placed from 
Solar System observations of planetary precessions and 
from laboratory experiments. While the LISA con- 
straints would beat those from Solar System observations 
(which presently give |aa| ;$ 1-2 x 10 18 m 2 ), considerably 
stronger constraints have already been placed from fifth- 
force tests[[3 Using existing results from the Eot-Wash 
experiment, we can constrain \a,2\ ^ 2 x 10~ 9 m 2 . For 
this range of 02, we would not expect the propagating 
Ricci mode to be excited by astrophysical systems as the 
cutoff frequency is too high. But, even in the absence of 
excitation of the Ricci mode, gravitational radiation in 



The LISA constraint relies upon the assumption that the weak- 
field metric does describe the exterior of a BH; there is no such 
caveat on the Solar System constraint since the weak-field metric 
is undoubtedly applicable for the spacetime exterior to the Sun. 



/(i?)-gravity is still modified through the dependence of 
the transverse polarizations on the Ricci scalar. 

Although the constraints from astrophysical observa- 
tions will be much weaker than this laboratory bound, 
they are still of interest since they probe gravity at a 
different scale and in a different environment. It is possi- 
ble that /(ii)-gravity is not universal, that it is different 
in different regions of space or at different energy scales. 
We could regard the f(R) model as an approximate ef- 
fective theory, and argue that the range of validity of a 
particular parameterization is limited to a specific scale. 
For example, we could imagine that the effective theory 
in the vicinity of a massive BH, where the curvature is 
large, is different from the appropriate effective theory in 
the Solar System, where curvature is small; or /(i?) could 
evolve with cosmological epoch so that it varies with red- 
shift. The limit on ai from gravitational- wave observa- 
tions will depend upon the BH mass, orbital radius and 
observation time, but it is clear that if the laboratory 
bound is indeed universal there should be no detectable 
deviation: observation of a deviation would thus prove 
not only that GR failed, but that the effective 02 varied 
with environment. 

One method of obtaining a variation is via the 
chameleon mechanism, where /( -R)-gravity is modified in 
the presence of matter flQ5WI07| . In metric f(R)- gravity 
this is a nonlinear effect arising from a large departure 
of the Ricci scalar from its background value [40( . The 
mass of the effective scalar degree of freedom then de- 
pends upon the density of its environment [57l Il08l ] . In 
a region of high matter density, such as the Earth, the 
deviations from standard gravity would be exponentially 
suppressed due to a large effective T; while on cosmo- 
logical scales, where the density is low, t he scalar would 
have a small T , perhaps of the order Ho / c |l05l . Il06| . The 
chameleon mechanism allows /(-R)-gravity to pass labo- 
ratory, or Solar System, tests while remaining of interest 
for cosmology. In the context of gravitational radiation, 
this would mean that the Ricci sc alar m ode could freely 
propagate on cosmological scales |l09l |. Unfortunately, 
since the chameleon mechanism suppresses the effects of 
f(R) in the presence of matter, this mode would have to 
be excited by something other than the acceleration of 
matter. Additionally since electromagnetic radiation has 
a traceless energy-momentum tensor it cannot excite the 
Ricci mode03 To be able to detect the Ricci mode we 
must observe it well away from any matter, which would 
cause it to become evanescent: a space-borne detector 
such as LISA could be our only hope. 



The standard transverse polarizations of gravitational radiation 
have an energy-momentum tensor that aver ages to be traceless, 
although this may not be the case locally |lld |: the contribu- 
tion to the gravitational averaged energy-momentum tensor from 
a propagating Ricci mode does have a nonzero trace, sec (1861 . 
In any case it is doubtful that gravitational energy-momentum 
could act as a source for detectable radiation. 



17 



As the chameleon mechanism is inherently nonlinear, it 
is difficult to discuss in terms of our linearized framework. 
Treating f(R) as an effective theory, we could incorporate 
the effects of matter by taking the coefficients {a n } to be 
functions of the matter stress-energy tensor (or its trace). 
In this case, the results presented here would hold in 
the event that the coefficient a 2 is slowly varying, such 
that it may be treated as approximately constant in the 
region of interest. The linearized wave equations, (f28|) 
and retain the same form in the case of a variable 
ci2, the only alteration would be that a 2 R^ replaces R^ 
as subject of the Klein-Gordon equation. In particular, 
the conclusion that 7 = 1 is unaffected by the possibility 
of a variable 0,2- 

An interesting extension to the work presented here 
would be to consider the case when the constant term 



in the function f(R), o,q, is nonzero. We would then be 
able to study perturbations with respect to (anti-)de Sit- 
ter space. This is relevant because the current ACDM 
paradigm indicates that we live in a universe with a pos- 
itive cosmological constant [13, lllll ] . Such a study would 
naturally complement an investigation into the effects of 
background curvature on propagation. 
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